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Abstract 



The ground state of a classical two-dimensional (2D) system with a finite 
number of charges particles, trapped by two positive impurity charges local- 
ized at a distance (zq) from the 2D plane and separated from each other by a 
distance x p are obtained. The impurities are allowed to carry more than one 
positive charge. This classical system can form a 2D-like classical molecule 
that exhibits structural transitions and spontaneous symmetry breaking as 
function of the separation between the positive charges before it transforms 
into two 2D-like classical atoms. We also observe structural transitions as a 
function of the dielectric constant of the substrate which supports the charged 
particles, in addition to broken symmetry states and unbinding of particles. 
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I. INTRODUCTION 



Classical charged particles confined in a two-dimensional (2D) layer has been a subject of 
interest in the last few years. |l]-[7| It is well-known that classical 2D electrons form a Wigner 
crystal. [l|,|J On the other hand, confined charged particles in 2D form 2D-like classical atoms 
which can exhibit very different melting and structural phase transitions. [|IIJ Recently 
these properties were studied for a set of charged particles confined by a parabolic and 
hard-wall well. PHT2[ They order in a ring structure and a Mendeleev table for these atomic 
like structures was constructed. An essentially different confinement is the one produced 
by a Coulomb potential. 2D charges which are confined by this potential also form 2D-like 
classical atomic structures but they exhibit very different structural transitions as compared 
to the one confined by parabolic potentials. |13| 

In the present paper we extend our previous work for single classical atoms ||,[13| and ver- 



tically coupled atoms |14| to laterally coupled atoms which form molecular structures. The 



system we consider here is composed of classical charged point particles interacting through 
a repulsive Coulomb potential and held together by two remote impurities a distance d away 
from the 2D plane where the charged particles can move. The ground state of the system is 
calculated by using the Monte Carlo simulation technique and by minimization techniques. 
We show that the ground state of this 2D-like classical molecule (2DCM) presents structural 
phase transitions and if the impurities are separated further apart it is broken into two 2D- 
like classical atoms (2DCA). Spontaneous broken symmetry and evaporation is observed, 
which depend strongly on the effective impurity charges and the distance between them. 

II. THEORETICAL MODEL 

Our model consists of a finite number of particles N, each of them with negative charge 
e = — | e |, which move in a 2D plane. Two fixed impurities, each of them with charge +Ze, 
are localized a distance z Q from the 2D plane and are separated by a distance Xp along 
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the x-axis and are placed inside a medium with dielectric constant e. All particles interact 
through a Coulomb potential (see Fig. 1). The energy for such a classical system is given 
by 

-Ze 2 N 2 1 e 2 N 1 

E = —J2E - - +~ E - - > (!) 

e i=i j=i | -Rj — -Dj | £o j>j=i | Ri — Rj | 

where Ri is the position of the i th negative particle at the (X, Y) plane and Dj = D x e x — Z Q e z 
(D x = Xp/2) is the position of the j th positive charge, i.e. the impurity. To keep the problem 
as simple as possible we have neglected in Eq. flU) the effects due to the image charges. 
Following our previous work [[H|] we will vary the dielectric constant e through which it is 
possible to change the effective confinement potential in a continuous way. For convenience, 
we express the energy in units of E Q = e 2 /e zo where e Q is the dielectric constant of vacuum 
and all the distances are now in units of z . This reduces Eq. ([!]) to 

_y* N 2 1 N -i 

c i=i j=i I rj — dj | j>i=i I '« 

where = E/E , Z* = Z/e*, x p = X P J z , dj = Dj/z , fj = Rj/z , e* = e/e . In order to 
obtain the minimum energy configuration we perform a numerical minimization of Eq. 
with respect to the 2N variables (coordinates (x,y) of the negative particles). 

III. RESULTS 

To perform the numerical simulation, we start from a random distribution of negative 
particles and obtain the minimum energy configuration by using standard minimization 
techniques. Next we change, e.g. x p , slightly and take as the initial configuration the last 
minimum obtained by our numerical minimization. To have an independent check we also use 
the Monte Carlo simulation technique starting from several different initial configurations. 
Different from the results obtained with the 2D classical atom |H| we found that bounded 



configurations , i.e., < oo V i, can exist if the total effective charge Zt = 2Z* satisfies 
the condition Zt < N — I for N > 1. In fact, we observed that the condition for bounded 



configurations tends to behave as a single 2-D classical atom, if the distance between the 
two positive charges goes to zero. To analyze the molecule-atom transition we fix the value 
of the dielectric constant e* = 1, we consider a system of six electrons and put the two 
positive charges at the points (d x , 0, —1) and (— d x , 0, —1), each with a charge Z = 3. In Fig. 
2 we show the distance of the six negative particles to the point (x, y, z) = (0, 0, 0) in the 
minimum energy configuration as a function of the distance between the positive charges 
x p = 2d x . We clearly observe four distinct regions. In the first region (I), x p is very small and 
the negative particles are arranged in two equilateral triangles, inside each other (see Fig. 
3(a)). The two positive charges (their positions are indicated by the crosses in Fig. 3) are so 
close to each other that they behave as one central charge. This configuration is very different 
from the one which was found in the case of a parabolic confinement potential || where one 
particle would be at (x,y) = (0,0) and the other 5 in a ring around it. At x p = 0.0884 the 
system undergoes a small structural transition to region (II) (see Fig. 3(b)) which is of first 
order. Note that increasing x p from x v = acts as a symmetry breaking field which breaks 
the rotational symmetry of the configuration depicted in Fig. 3(a). With increasing x p , 
the ground state configuration changes continuously up to x p = 1.0074 where it undergoes 
a structural transition and a new configuration appears in region (III). This transition is 
characterized by the fact that the particles change their positions abruptly (compare solid 
(II) and open (III) symbols in Fig. 3(c)). With further increase of x p the system exhibits 
a new structural transition to region (IV) at x p = 1.8792 where the configuration clearly 
exhibits now inversion symmetry with respect to the center of the coordinate system. In 
region (IV) the negative particles tend to stay in two triangles, each one close to one of the 
positive charges. Up to this region, the system behaves as a classical molecule. Namely, all 
distances between the negative particles are related to the positions of both positive charges. 
Finally, for a value x p ~ 3.0 we observe more explicitly that the negative particles split in 
two sets, each one correlated to one positive charge. No structural transition is associated 
with this behavior but, up to this point, the system behaves as two independent 2D classical 
atoms. This can be seen by the fact that, up to x p ~ 3.0 the system exhibits rotational 
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symmetry along the z-axis for each of the two positive charges. Different from the structural 
transition from regions (I) up to (IV), the one from 2D classical molecule to 2D classical 
atom is not characterized by a structural transition, i.e. it is continuous. In order to visualize 
these stable states, in Fig. 3 we show typical configurations of the system in each region 
and in Fig. 4 we plot the distance of three negative particles related to each positive charge. 
As shown in Fig. 4, when x p > 3.0 the system behaves as two independent atoms, with the 
negative particles in two equilater triangles and the positive charges in the center of each 
one, with energy E sing i e — > —4.899. In order to show the order of the structural transitions 
we calculated the first of the energy as a function of x p , by numerically differentiating the 
energy curve. The results for the first derivative are presented in Fig. 5. As can be seen 
the first derivative is discontinuous at the structural transitions (I — > II, II — > III and III — > 
IV). It is continuous throughout the region where the classical 2D molecule transits to two 
2D classical atoms. 

Different structural transitions are observed if we analyze the behavior of the minimum 
energy configuration as a function of the dielectric constant (e*). To this, we consider two 
positive charges (each one with Z — 3) fixed at the points (1,0,-1) and (—1,0,-1) and 
six negative particles which are free to move in the xy-plane. For e* = 1 we are in region 
IV of the above. In Fig. 6 we show the distance of those six negative particles to the 
point (x,y,z) = (0,0,0) for the minimum energy configuration as a function of e*. Eleven 
different regions are observed, of which we plot in Fig. 6 only the first ten of them. In 
region (I), 1.0 < e* < 1.054, the minimum energy configuration corresponds to the negative 
particles arranged in two triangles each one close to each positive charge (see Fig. 7(a)). 
When e* is increased a structural transition takes place at e* = 1.117. Up to this value of e*, 
region (III), the configuration changes continuously and for e* — > 1.225 one of the negative 
particles becomes weakly bound and moves to infinity, i.e. it evaporates. After this value, 
region (IV), the minimum energy configuration changes continuously and again presents a 
structural transition, at e* = 1.314, to a small region (V). This region goes up to e* = 1.325, 
where a new structural transition takes place. Now, in the new region, region (VI), another 
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negative particle starts to become unbound when the value e* = 1.492 is reached. The 
minimum energy configuration is stable in the new region (VII) until the value e* = 1.909 
and again changes its structure. After this point, again another particle starts to move 
away, region (VIII), and goes to infinity at e* = 2.112. Up to this point, region (IX), only 
three particles remain bound until e* = 3.000. In region(X) 16.2 > e* > 3.000 only two 
particles remain, with each negative particle sitting almost on top of each of the positive 
charges in a symmetrical arrangement. Finally, when e* > 16.2 only one particle remains 
bounded. Unexpectedly, this negative particle stays in an asymmetric position, in this case 
at a distance r\ = 0.77. This fact can be confirmed by taking the minimum of the potential 
energy. In this case it presents a minimum at r\ = 0.77 and a local maximum at r\ = 0, 
which corresponds to a symmetrical arrangement. The existence of region (X) shows that 
a bounded configuration can exist in a region where the condition < N — 1 for N > 1 is 
not satisfied. 

IV. CONCLUSIONS 

The molecule-atom transition and structural transitions were studied in a 2D model sys- 
tem of negative particles confined by a Coulomb potential of two positive impurity charges. 
This system presents a large variety of structural transitions as a function of its parameters, 
the distance between the impurity charges and the dielectric constant. As observed for a 2D 
like atom the 2D like molecule exhibits highly asymmetric ground state configurations which 
are related to broken symmetry states. However, the condition Z T < N — 1 for the existence 
of a bounded configuration is only satisfied if the positive particles are very close. Finally, 
our results show that a finite 2D system confined by a nonlinear potential presents a much 
larger wealth of structural transitions and new states as compared to the ones confined by 
parabolic potentials. 
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FIGURES 



FIG. 1. Schematic representation of the system considered in this paper. 

FIG. 2. Distance ri of the six negative particles to the point (x, y, z) = (0, 0, 0) for the minimum 
energy configuration as a function of the distance between the two positive particles x p . Structural 
transitions are indicated by the vertical dotted lines. 

FIG. 3. Typical configurations (and transitions) of the system in the different regions corre- 
sponding to Fig. 2. 

FIG. 4. Distance of three negative particles related to each positive charge. Structural transi- 
tions are indicated by the vertical dotted lines. 

FIG. 5. The first derivative of the energy as a function of x p . The inset is an enlargement of 
the small x p region of the discontinuity of the first derivative. Structural transitions are indicated 
by the vertical dotted lines. 

FIG. 6. Distance of six negative particles to the point (x, y, z) = (0, 0, 0) in the minimum energy 
configuration as a function of the logarithm of the dielectric constant e*. The dashed lines indicate 
structural transitions while the vertical dotted lines correspond to the evaporation of a particle. 

FIG. 7. Typical stable configurations of the ground state in different regions of e*. 
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This figure "fig 1. PNG" is available in "PNG" format from: 



http://arXiv.org/ps/cond-mat/0203137vl 




Fig. 3 
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